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ABSTRACT 

Recent observations suggest that helium became fully ionized around redshift z ~ 3. The He II 
optical depth derived from the Lyman-a forest decreases substantially from this period to z ~ 2; 
moreover, it fluctuates strongly near z ~ 3 and then evolves smoothly at lower redshifts. From these 
opacities, we compute, using a semi-analytic model, the evolution of the mean photoionization rate 
and the attenuation length for helium over the redshift range 2.0 < z < 3.2. This model includes an 
inhomogeneous metagalactic radiation background, which is expected during and after helium reion- 
ization. We find that assuming a uniform background underestimates the required photoionization 
rate by up to a factor ~ 2. When averaged over the (few) available lines of sight, the effective optical 
depth exhibits a discontinuity near z w 2.8, but the measurement uncertainties are sizable. This 
feature translates into a jump in the photoionization rate and, provided the quasar emissivity evolves 
smoothly, in the effective attenuation length, perhaps signaling the helium reionization era. We then 
compute the evolution of the effective optical depth for a variety of simple helium reionization models, 
in which the measured quasar luminosity function and the attenuation length, as well as the evolving 
He III fraction, are inputs. A model with reionization ending around redshift z m 2.7 is most consistent 
with the data, although the constraints are not strong thanks to the sparseness of the data. 
Subject headings: intergalactic medium - diffuse radiation - quasars: absorption lines 



1. INTRODUCTION 

In the standard reionization history of the Universe, 
the metagalactic ionizing background evolved relatively 
slowly except for the reionization of hydrogen (z > 6) 
and of helium (fully- ionized at z ~ 3). During these 
phase transitions, the intergalactic medium (IGM) be- 
came mostly transparent to the relevant ionizing pho- 
tons, allowing the high-energy radiation field to grow 
rapidly. Since He II has an ionization potential of 
54.4 eV, bright quasars with hard spectra are required 
for its ionization. As such, the distribution and intrinsic 
character of quasars, in addition to the properties of the 
IGM, determine the radiation background at these high 
energies. 

These quasars are quite rare, implying a strongly fluc- 
tuating background even after reionization dFardal et al.l 
119981 : iBolton et alJl2006t lMeiksinll2007t lFurlanettoll2008f ) ■ 
Direct evidence for these source-induced variations has 
been seen in the "transverse proximity effect" of the 
hardness ratio through comparisons of the H I and 
He II L yman-a (Lya) forests with surveys for nearby 
Jakobsen et al.ll2003HWorseck fc Wisotzkil200a 
These variations are exag- 
attenuation from the IGM 



qua sars 
iWorseck et al.l 120071) . 

aerated by the stron; 



gerated by trie strong attenuatio n trom tne HjM 
dHaardt fc Madaul 119961 : iFaucher-Giguere et all l2008aL 
iFurlanetto fc Obi l2008al ) . Furthermore, radiative trans- 
fer through the clumpy IGM can induce additional 
fluctu ations (jMaselli fc Ferraral [20051 : iTittlev fc Meiksinl 
2007). During reionization, fluctuations in the back- 
ground are even greater, because some regions receive 
strong ionizing radiation while others remain singly- 
ionized with no local illumination. 

Recent observations indicate that helium reionization 
occurs at z ~ 3. The strongest evidence comes from far- 
ultraviolet spectra of the He II Lya forest along the lines 
of sight to bright quasars at z ~ 3. These observations 



of the He II Lya transition (A rcs t = 304 A) are difficult, 
because bright quasars with sufficient far-UV flux and no 
intervening Lyman-limit systems are required. To date, 
six such lines of sig ht have yielded opac it y measurements: 
PKS 1935-692 dTytler et al.l Il9p lAnderson et al " 



1999)), HS 1700+64 (iDavidsen et al.lll996tlFechner et al 
2006) HE 2347-4342 dReimers et all 119971; iKriss eta! 
boOlUSmette et al.l 12002b IShull et al . 2004: Z heng et al 
2004bh . SP SS .12346-0016 (IZheng et all l2004ai l2008h 



Q0302-003 (Dakobsen et all Il994t IHoe 
iHeap et al.l 120001: Uakobsen et al.H 2003l) 



iHogan et alj 119971: 
„ and HS 1157- 

3143 (jReimers et al.ll2005l ). The effective helium optical 
depth from these studies decreases rapidly at z rj 2.8, 
then declines slowly to lower redshifts. The opacities at 
higher redshifts exhibit a patchy structure with alter- 
nating high and low absorption, which may indicate an 
inhomogeneous radiation bac kground. More prom ising 
sightlines have been detected (jSvphers et alj l2009h . and 
the Cosmic Origins Spectrograph on the Hubble Space 
Telescope should add to the current pool of data. 

Several indirect methods attempt to probe the im- 
pact of helium reionization on the the IGM. One ex- 
pected effect of heliu m reionization is an increase in 
the I G M temperature dHui fc Gnedinlll997tlGleser et all 
s Or]"" 



20051; IFurlanetto 
Schave et all (|2000f ) 

crease at z ~ 3.3 



h008bl : iMcQuinn et all 12009) 
detected a sudden temperature in- 
by examining the th ermal broad- 
ening of H I Lya forest lines (see also ISchave et alj 
119991 : iTheuns et all l2002bl ). Around the same time, 
the IGM temperature-density relation appears to be- 
come nearly isother mal, another i ndica ti on of recent 
helium reionization (|Schave et alj 120001: iRicotti et a! 



2000). However, not all studies agree (jMcDonald et al 
|2001| ). and temperature measurements via the Lya for- 
est flux powe r spectrum show no evidence for any sud- 
den change (jZaldarriaga et all 120011 : IViel et all 120041 : 



2 



iMcDonald et al.ll2006l ). Furthermore, this temperature 
increase should decrease the reco mbination rate of hy - 
drogen, decreasing the H I opacity (|Theuns et al.ll2002af ). 
Three studies with differing methods have measured 
a narrow dip at z ~ 3.2 in the H I effective opti- 
cal de pth dBernardi et alJ 120031 : iFaucher-Giguere et al.1 
I2008U iDall' Aglio et al II2008T). While initia lly attributed 
to helium reionization ( Theuns et al.lr2002af ). more recent 
studies find that reproducing this feature with helium 
reionization is ex t remely difficult (iBolton et al.l 120091 : 
iBolton et alj|2009t iMcQuinn et al.ll2009D . 

The (average) metagalactic ionizing background 
should also harden as helium is reionized, because the 
IGM would become increasingly transparent to high- 
energy photons. Observations of the He II/H I ratio 
are qualitatively consistent with reio nization occurrin g 
at z ~ 3 in at l east o ne line of sight (|Heap et al.ll2000D . 
ISongailal (|1998l l2005f) found a break in the ratio of C IV 
to Si IV at z ~ 3. Modeling of the ionizing background 
from optically thin and optically thick met al line systems 
also shows a significant hardening at z ~ 3 (|Vladilo et al.l 
120031: lAgafonova et~ai] l2005l I2007D . but other data of 
comparable qual i ty show no evidence for rapid evolution 
(|Kim et alj|2002t lAguirre et al.ll2004h . However, this ap- 
proach is made more difficult by the large fluctuations 
in the He II/H I ratio even after helium reionization is 
complete (|Shull et al.ll2004l ). 

In this paper, we focus on interpreting the He II Lya 
forest and the significance of the jump in the opacity at 
z « 2.8. After averaging the effective optical depth over 
all sightlines, we calculate the expected photoionization 
rate given some simple assumptions. In particular, we 
investigate the impact of a fluctuating radiation back- 
ground, comparing it to the common uniform assump- 
tion. We interpret our results in terms of an evolving at- 
tenuation length for helium-ionizing photons Rq as well 
as state-of-the-art models of inhomogeneous reionization. 

We use a semi-analytic model, outlined in to in- 
fer the helium photoionization rate from the He II Lya 
forest. The helium opacity measurements in the redshift 
range 2.0 < z < 3.2, which serve as the foundation for 
our calculations, are compiled from the literature in Ej3j 
First, we assume a post-reionization universe over the 
entire redshift span. In this regime, we find the pho- 
toionization rate and attenuation length in given the 
average measured opacity. Motivated by these results, 
we examine some fiducial reionization histories in Sj5j We 
conclude in Sj6] 

In our calculations, we assume a cosmology with Q rn = 
0.26, Cl A = 0.74, n b = 0.044, H = h(W0 km s" 1 Mpc" 1 ) 
(with h — 0.74), n — 0.95, and < t s = 0.8, consistent 
with the most recent measurements (|Dunklev et al.ll2009l : 
iKomatsu et~aT]|2009D . 

2. METHOD 

The helium Lya forest observed in the spectra of 
quasars originates from singly-ionized helium gas in the 
IGM. Quantitative measurement of this absorption is 
typically quoted as the transmitted flux ratio F, defined 
as the ratio of observed and intrinsic fluxes, or the related 
effective optical depth 

r eff = -ln(F). (1) 
We have F = e^ T <= tf = (e~ T ) > e~ (r \ thus r cff < (r). 
From the current opacity measurements, we wish to infer 



the He II photoionization rate. This connection depends 
on the details of the IGM, including the temperature, 
density distribution, and ionized helium fraction. 

2.1. Fluctuating Gunn-Peterson approximation 

The Gunn-Peterson |l965j) optical depth for He II Lya 
photons is 



TGP = faKH 1 (z)n H eII- 

m e c 



(2) 



Here, the oscillator strength f a — 0.416, X a — 304 A, and 
TiHeii is the density of singly- ionized helium in the IGM. 
For simplicity, we approximate the Hubble constant as 

Since the Lya forest probes 
the low-density, ionized IGM, most of the hydrogen (mass 
fraction X = 0.76) and helium (Y — 0.24) are in the form 
of H II and He III, respectively, after reionization. Under 
these assumptions, photoionization equilibrium requires 



rn H eII = "HeleQA, 



(3) 



where T is the He II photoionization rate 
and the case-A recombination coefficient is 
a A = 3.54 x 1Q- 12 (T/10 4 K )~ - 7 cm 3 s" 1 accord- 
ing to lStorev fc Humme r (1995). For a clumpy universe, 
most photons emitted by recombinations are produced 
and subsequently reabsorbed in dense, mostly neutral 
systems. These ionizing photons, therefore, cannot 
escape to the low dens ity regions of interes t for the 
forest, so we use case-A (|Miralda-Escudell2003f) . 

The Lya forest, and therefore the optical depth, trace 
the local overdensity A of the IGM, where A = p/p 
and p is the mean mass density. Since n e cx nuo in a 
highly ionized IGM, equation ([5]) implies that nndi °c 
n\ a cx A 2 . The optical depth is proportional to nneii (see 
cq. [2]), which introduces a A 2 factor. Additionally, the 
temperature of the IGM, which affects the recombination 
rate, is typica lly described by a p ower law of the form 
T = ToA't- 1 (jHui fc Gnedid 11991 . where To and 7 are 
taken as constants. 1 Including the above equations and 
cosmological factors, 

"■ 7 12 \ 2 fn m h 2 y 1/2 

7 V 0.142/ 



TGP ; 



19 



T 
10 4 K 



0.0241, 



1 + 2 



9/2 



A 



2-0.7(7-1) 



(4) 



where T = 10" 14 r_i 4 s _1 . 

The fluctua ting Gunn-Peterson approximation 
(FGPA) (e.g., IWeinberg etldl I1999I ). equation ©, 
relates the effective optical depth, or the continuum 
normalized flux, to the local overdensity A and the 
photoionization rate T_i4. This approximation shows 
the relationship between the opacity and the IGM, but 
it ignores the effects of peculiar velocities and thermal 
broadening on the Lya lines. In practice, when compar- 
ing to observations an overall proportionality constant, 
k, is introduced to the right hand side of equation (01 
to compensate for these factors, as described in <j2.3l 
This normalization also incorporates the uncertainties 
in To, fl m , and h. 

1 In actuality, during and after helium reionization both To and 
7 likely become redshift- and density-dependent. We will consider 
such effects in El and fl4"3l 
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2.2. A semi-analytic model for Lya absorption 

To fully describe helium reionization and compute 
the detailed features of the Lya forest, complex hy- 
drodynamical simulations of the IGM, including radia- 
tive transfer effects and an inhomog eneous background , 
are required. Recent simulations (ISokasian et al.l l2002t 
iPaschos et all l2007t iMcQuinn et all I2009D have made 
great advances to incorporate the relevant physics and to 
increase in scale. The simulations remain computation- 
ally intensive, and they cannot simultaneously resolve 
the ^100 Mpc scales required to adequately study in- 
homogeneous helium reionization and the much smaller 
scales required to self-consistently study the Lya forest, 
necessitating some sort of semi-analytic prescription to 
describe baryonic matter on small scales. On the other 
hand, the semi-analytic approach taken here, including 
fluctuations in the ionizing background, should broadly 
reproduce the observed optical depth, especially consid- 
ering the uncertainties in the measurements and the lim- 
ited availability of suitable quasar lines of sight. To 
outline, the model has four basic inputs: the IGM den- 
sity distribution p(A), the temperature-density relation 
T(A), the radiation background distribution f(J), and 
t he mean helium ionized fractio n xhcIH- 

iMiralda-Escude et al.l 1)20001 ) suggest the volume- 
weighted density distribution function 

(A- 2 /3 _ C ) 21 



p(A) = AA' 13 exp 



(5) 



2 (2<y3r 

where 5q — 7.61/(1 + z) and j3 (for a few redshifts) are 
given in Table 1 of their paper. Intermediate (3 values 
were found using polynomial interpolation. The remain- 
ing constants, A and Co, were calculated by normaliz- 
ing the total volume and mass to unity at each redshift. 
The distribution matches cosmological simulations rea- 
sonably well for the redshifts of interest, 2 i.e. z = 2 — 4. 
Although this form does not incorporate all the physics of 
reionization and was not generated with the current cos- 
mological parameter values, the overall behavior should 
be sufficient for the purposes of our model. 

A current topic of discussion is the thermal evo- 
lution of the IGM during and after helium reioniza- 
tion. The temperature-density relation should vary 
as a function of redshift and density; notably, helium 
reionization should i ncrease the overall temperature by 
a fac t or of a few dHui fe Gnedinl 119971: iGleser et al.l 
I2005t iFurlanetto fc Ohl l2008bt IMcQuinn et all 12009). 
which may have been observed (jSchave et al.l 120001 : 
iRicotti et al.ll200 0). As mentioned in §2.1[ the tempera- 
ture is assumed to follow a power law T = XoA 7-1 in the 
FGPA. Generally, T Q ~ 1 - 2 x 10 4 K and 1 < 7 < 1.6 
should broadly descr ibe the post-reionization IGM (e.g., 
iHui fe Gnedinl [T997f ). but the exact values are a matter 
of debate. Unless otherwise noted, we use To = 2 x 10 4 K 
and 7=1 throughout our calculations. The isothermal 
assumption also suppresses the temperature, and there- 
fore density, dependence of the recombination coefficient 
a a- 

A uniform radiation background has been a com- 
mon assumption in previous studies, but the sources 



(quasars) for these photons are rare and bright. There- 
fore, random variations in the quasar distribution cre- 
ate substantial v ariations in the high-frequency radia- 
tion background (lMeiksirJ l2007t IFurlanetto h Ohll2008al : 
IMcQuinn et al.ll2009h . Furthermore, the 1/r 2 intensity 
profiles of these sources induce strong small-scale fluctu- 
ations (|Furlanettol 12009) , which may in turn significantly 
affect the overall optical depth of the Lya forest. For the 
probability distribution /(J) of the angle- averaged spe- 
cific intensity of the ra diation backgroun d </, we follow 
the model presented in iFurlanettol ([2008) . In the post- 
reionization limit, the probability distribution can be 
computed exactly for a given quasar luminosity function 
and attenuation length, assuming that the sources are 
randomly distributed (following Poisson statistics). This 
distributio n can be derived either via Markov's method 
Zu^lll99l) or via the method of characteristic functions 



Meiksin fe White! 120031. 



2 More rece n t sim ulations by Pawlik et al. ( 2009;) and 
I Bolton fc Becker (2009) basically agree with the above p(A) for 
low densities, which the Lya forest primarily probes. 



During reionization, the local He III bubble radius, 
i.e. the horizon within which ionizing sources are visible, 
varies across the IGM and so becomes another impor- 
tant parameter. Due to the rarity of sources, with typi- 
cally only a few visible per He III region, a Monte Carlo 
treatment best serves this regime (|Furlanettdl2"008l ). For 
a given bubble of a specified size, we randomly choose 
the number of quasars inside each bubble according to 
a Poisson distribution. Each quasar is then randomly 
assigned a location within the bubble as well as a lu- 
minosity (via the measured luminosity function). Next, 
we sum the specific intensity from each quasar. After 
10 6 Monte Carlo trials, this procedure provides /(J) for 
a given bubble size; the solution converges to the post- 
reionization scenario for an infinite bubble radius. The 
final ingredient is the size distribution of discrete ion- 
ized bub bles, which is found using the excursion set ap- 
proach of lFurlanetto fc Ohl (|2008al) (based on the hydro - 
gen reionization equivalent from IFurlanetto et al.ll2004T ). 
After integrating over all possible bubble sizes, /(J) de- 
pends on the parameters z, Rq, and ShcIII, in addition 
to the specified luminosity function. 

In the following calculations, we scale J (which is eval- 
uated at a single frequency) to T (which integrates over 
all frequencies) simply using the ratio of the mean pho- 
toionization rate to the mean radiation background. This 
is not strictly correct, because higher frequency pho- 
tons have larger attenuation lengths and so more uniform 
backgrounds; however, it is a reasonable prescription be- 
cause the ionization cross section falls rapidly with pho- 
ton frequency. However, it does mean that we ignore the 
large range in spectral indices of the ionizing sources (see 
below), which modulate the shape of the local ionizing 
background and lead to an additional source of fluctu- 
ations in r relative to J that we do not model. The 
largest problem occurs during helium reionization, when 
the highest energy photons can travel between He III re- 
gions (a process we ignore); however, they have small ion- 
ization cross sections and so do not significantly change 
our results, exce pt very near the end of that process 
(|Furlanettoll2008l ). 

For the majority of the paper, we consider the IGM 
to be fully-ionized, i.e. ShcIII = 1-0. In this post- 
reionization regime, Rq, the attenuation length of the 
ionizing photons, determines the shape of /(r), given 
the redshift and other model assumptions. We show 
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Fig. 1. — Distribution of the photoionization rate T relative to its mean value in a fully-ionized IGM, (r), at redshift z = 2.5. Left panel: 
The curves assume Ro = 5, 10, 50, and 100 Mpc, from widest to narrowest, in a post-reionization universe. Right panel: The curves take 
^Helll = 0.3,0.5,0.75,0.9, and 1.0, from lowest to highest at peak, for Ro = 35 Mpc. 



several example distributions in the left panel Figure [TJ 
From widest to narrowest, the curves correspond to 
R = 5, 10,50, and 100 Mpc (z = 2.5). Smaller attenu- 
ation lengths yield a greater spread in T. Qualitatively, 
more sources contribute to ionizing a given patch of the 
IGM for higher Ro, making the peak photoionization rate 
more likely, i.e. the curve is narrower. A uniform back- 
ground corresponds to Ro — > oc. Although low T values 
are more likely for low Ro, the high-r tail is nearly inde- 
pendent of Ro- This is because a large T occurs within 
the "proximity zone" of a single quasar, making it rel- 
atively independent of contributions from much larger 
scales (unless Ro is much smaller than the proximity zone 
itself) . During reionization, as in Ej5j the ionized fraction 
and the mean free path affect the distribution function, 
as shown in the right panel of Figure [TJ Here, we take 
z = 2.5, Ro = 35 Mpc, and XHein = 0.3, 0.5, 0.75, 0.9, and 
1.0 (from lowest to highest at the peak). A broader distri- 
bution of photoionization rates is expected, because the 
large spread in He III bubble sizes restricts the source 
horizon inhomogencously across the Universe. 

To estimate r g (or F), we integrate over all densities 
and photoionization rates: 



rate is 



47T 



—a v dv, 
hv 



(7) 



where cj„ = 1.91 x 10 _18 (^/t , Heii) -3 cm 2 is the photoion- 
ization cross section for He II and i^Hcii is the photon fre- 
quency needed to fully ionize helium. For the radiation 
background at frequency v, J v , we as sume a simplified 
form, the absorption limited case in iMeiksin fc White! 
(|200l : 

J v = ^-e y {z)Ro{z), (8) 

47T 

where e„ is the quasar emissivity and Ro is the attenua- 
tion length. 

We begin with the B-band emiss ivity eg , derived from 
the quasar luminosity function in [Hopkins "eTall ((2007I . 
hereafter HRH07). To convert this to the extreme-UV 
(EUV) frequencies of interest, we follow a broken pq wer- 
law spectral energy distribution (|Madau et al.| [l999): 



L{y) oc < v 



-0.3 
-0.8 



2500 < A < 4600 A 
1050 < A < 2500 A 
A < 1050 A. 



(9) 



F 



Jo Jo 



(6) 



This integral is valid if T and A are uncorrelated. Since 
relatively rare quasars ionize He II, random fluctua- 
tions in the number of sources (as opposed to their 
spatial clustering) dom inate the ionization morphology 
(|McQuinn et al.lT2009t ). justifying our assumption of an 
uncorrelated IGM density and the photoionization rate. 

2.3. The UV background from quasars 

Since quasars ionize the He II in the IGM, the UV 
metagalactic background can, in principle, be calcu- 
lated directly from distribution and intrinsic properties 
of quasars. Currently, these details, i.e. the quasar lu- 
minosity function and attenuation length, are uncertain. 
The following method for estimating T is used as a ref- 
erence for the semi-analytic model. The He II ionization 



The EUV spectral index a is a source of debate and is 
not the same for all quasars. iTelfer et alj (|2002[ ) find a 
wide range of values for individual quasars, e.g. a = -0.56 
for HE 2347-4342 and 5.29 for TON 34. Most quasars 
lie closer to the mean, but HE 2347-4342 is a He II Lya 
line of sight in iJS] Unless otherwise noted, w e use the 
mean value (a) 1 6 fromlTelfer et all (|2002f) , ignoring 
the variations in a. iZheng et all (119971 ) found (a) w 1.8, 
which serves as a comparison. The uncertainty in the 
spectral index a affects the amplitude, not the shape, of 
the r curve derived from the QLF (see Fig. [2]). Since 
our semi-analytic calculations are normalized to a single 
point on the T curve, this uncertainty translates into an 
amplitude shift in our results, i.e. changes k. 

The next ingredient for J„ is the attenuation length of 
helium- ionizing photons Ro- As described above, Ro de- 
pends on the photon frequency; for example, high-energy 
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Fig. 2. — Evolution of the mean photoionization rate (in units of 
10~ 14 s^ 1 ) with redshift. The solid curve represents the inferred 
ionization rate from the HRH07 quasar lumi nosity function with 
a = 1.6. The dotted (dashed) curve follows Madau et all (IT999T) 
with extreme-UV spectral ind ex a = 1.6 (1.8). All take the at- 
tenuation length from eq. (ID) . The points (with a slight redshift 
offset) show the reconstructed photoionization rate given the mea- 
sured effective optical depth, assuming a uniform (triangles) and 
fluctuating (squares) radiation background. The results are fixed 
to the solid curve at z = 2.45. 

photons can propagate larger distances. For simplicity, 
we use a single frequency-averaged attenuation length 
and focus only on the redshift evolution (in any case, the 
absolute amplitude can be subsumed into our normal- 
ization factor k below). For concreteness, w e apply the 
comoving form found in lBolton et alj (|2006f ): 

R = 30 M P C : (10) 

which assumes the number of Lyman limit sys- 
tems per unit redshift is p roportional to (1 + z) 1 - 5 
(jStorrie-Lombardi et al.l I1994D and uses the normaliza - 
tion based on the model o f lMiralda-Escude et al.l (|2000D . 
An alternate approach, probably more appropriate dur- 
ing reionization itself, is to estimate the attenuation 
length around individual quasars, as in lFurlanetto fc Obi 
(2008a|). This method gives a similar value at z = 3 
but with a slower redshift evolution, which would only 
strengthen our conclusions. This empirical equation for 
the mean free path (with the above quasar emissivity) 
provides the photoionization rate r , shown in Fi g ure [2] 
The photoionization rate from the iMadau et al.l (|1999f ) 
QLF (with a = 1.6 and 1.8) is also included in the fig- 
ure to illustrate the effect of uncertainties in the quasar 
properties. The spectral index has a greater effect on the 
amplitude than the differing QLFs. 

As noted in H2.ll comparing our FGPA treatment to 
the real Lya forest data requires an uncertain correction 
to the t — A relation in equationJU For this purpose, we 
assume the above emissivity and attenuation length to 
be accurate. Then, the semi-analytic model is adjusted 
so that the predicted T e s matches the measured value at 
a particular redshift. To do so, we insert a prefactor, 
k, to the right hand side of equation (JJ). This factor 
compensates for line blending and other detailed physics 



ignored by the FGPA, but it also includes any uncertain- 
ties in the underlying cosmological or IGM parameters 
(such as To). A suitable normalization redshift should 
be after reionization and have data from more than one 
line of sight (see the right panel of Fig. [3|). Throughout 
this work, we take z = 2.45 as our fiducial point. 

If k is a constant with redshift, our choice of refer- 
ence point mainly affects the overall amplitude of the 
photoionization rate or attenuation length, not the red- 
shift evolution. However, since k depends on the IGM 
properties, it may change with redshift, density, and/or 
temperature. For example, an increase in temperature 
broadens the widths of the absorption lines in the Lya 
forest, decreasing the importance of saturation but in- 
creasing the likelihood of line blending. Reionization, a 
drastic change to the IGM, should also affect k. Here, 
we take k to be independent of z. The precise value lies 
somewhere between 0.1 and 0.5, de pending on the spe- 
cific model. iFurlanetto fc Ohl (|2009l ) find k = 0.3 for the 
post-reionization hydrogen Lya effective optical depth, 
which is similar to our values. In any case, we emphasize 
that our method cannot be used to estimate the absolute 
value of the ionizing background - for which detailed sim- 
ulations are necessary - but we hope that it can address 
the redshift evolution of T. 

3. EVOLUTION OF THE HE II EFFECTIVE OPTICAL 
DEPTH 

Measurements of the He II effective optical depth are 
challenging. Suitable lines of sight require a bright 
quasar with sufficient far-UV flux and no intervening 
Lyman-limit systems. Currently, only five quasar spectra 
have provided He II opacity measurements appropriate 
for our analysis, displaye d in the left panel of Figure [3] 
with the averaged values. IZheng etTa l. (2004a) measure 
a lower limit on the optical depth at z ~ 3.5 from a sixth 
sightline, SDSS J2346-0016; we do not include any lower 
limits in any subsequent calculations, but we reference 
this measurement in Due to the limited scope of the 
data, the observed opacities may not be representative 
of the IGM as a whole. 

HE 2347-4342: This quasar (z cm = 2.885) is especially 
bright; ther efore, this lin e of sigh t has been extensively 
analyzed. IZheng et al.l (|2004bl ) completed the most 
comprehensive investigation, covering the r edshift range 
2.0 < z < 2.9, including Lya and Ly/3. IZheng et al.l 
(|2004bfl and lShull et all pOOl (2.0 < z < 2.9 and Lya 
only) utilized high-resolution spectra from the Far Ul- 
traviolet Spectroscopic Explorer (FUSE; R ~ 20,000) 
and the Very Large Telescope ( VLT; R ~ 45,000 ). An 
older and lower resolution study, Kriss et al. (2001), cov- 
ered 2.3 < z < 2.7 with FUSE. Below redshift z = 2.7, 
the effective helium optical depth evolves smoothly. At 
higher redshifts, the opacity exhibits a patchy structure 
with very low and very high absorption, often described 
respectively as voids and filaments in the literature. 

HS 1700+64: The Lya forest of this quasar (z cm = 
2.72) has been resolved with FUSE over t he redshift 
range 2.29 < z < 2.75 (jFechner et al.ll2006l ). An older 
study using the Hopkins Ultraviolet Telescope (HUT) 
is consistent wi t h the newer, higher resolution results 
(jDavidsen et al.l Il996h . The helium opacity evolves 
smoothly and exhibits no indication of reionization. 

Q0302-003: The spectrum of this quasar (z cm = 3.286) 
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Fig. 3. — Evolution of the He II effective optical depth based on the observations of the Lya forest for five quasar spectra: HE 2347-4342, 
HS 1700+64, Q0302-003, HS 1157+314, and PKS 1935-692. The squares are the opacity measurements averaged over redshift bins Az = 0.1 
with suggestive uncertainties. Left panel: Data and uncertainties as quoted in the literature are plotted. Right panel: The average values 
for each line of sight are displayed, elucidating the origin of the uncertainties in the average opacities used throughout the paper. The small 
redshift offsets within each bin are for illustrative purposes only. 



was observed with the Space Telescope Imaging Spec- 
trograph (STIS) aboard the Hubble Space Telescope 
(HST) at 1.8 A resolution (|Heap et alj|200"ch . The ef- 
fective optical depth generally increases with increas- 
ing redshift over 2.77 < z < 3.22, exclud ing a void at 
z ~ 3.05 due to a nearby ionizing source (iBaitlik et ail 
fl988tlZheng fc Davidsenill995l : lOiroux et aLlflQQl . The 
data were average d over redshift bins of A 2 ~ 0.1. 
iHogan et~a l. (1997) presented an analysis using the God- 
dard High Resolution Spectrograph (GHRS) which gen- 
erally agrees with the later study but quoted a noticeably 
lower He II optic al depth near z ~ 3.1 5. 

HS 1157+314: iReimers et all (|2005t l obtained low res- 
olution HST/STIS spectra of the He II Lya forest to- 
ward this quasar (z cm <~ 3). Over the redshift range 
(2.75 < z < 2.97) of the study a patchy structure, simi- 
lar to HE 2347-4342, is present. The given optical depth 
was averaged over a redshift bin of Az ~ 0.1. 

PKS 1935-692: The HST/STIS sp ectrum for this 
quasa r (z cm = 3.18) was analyzed by lAnderson et alJ 
(1999). Only one optical depth was quoted, but the spec- 
trum exhibited the usual fluctuations. 

To merge these data sets, we initially binned them in 
redshift intervals of 0.1, starting with z = 2.0. Each bin 
was assigned the median redshift value, e.g. z — 2.35 
for 2.3 < z < 2.4. To objectively combine the data, 
the transmission flux ratios for each data set were aver- 
aged (weighted by redshift coverage) in the redshift bins. 
Then, the values for each quasar were averaged, since 
multiple data sets may cover the same line of sight. Fi- 
nally, if more than one line of sight contributes to a bin, 
the fluxes are averaged once again. The process trans- 
lates the left panel to the right panel of Figure [3] The 
uncertainties along each line of sight in the right panel are 
simply the errors from each separate point in the litera- 
ture, added in quadrature, without regard to systematic 
errors; our averaged values then take errors equal to the 
range spanned by these separate lines of sight. For the 



remainder of the paper, the error bars result from assum- 
ing F ± dF for each redshift bin, where F are the squares 
and ±dF are the upper/lower error bars on the squares 
in the figure. The small number of well-studied lines of 
sight limits the amount of truly quantitative statements 
that can be made. 

4. RESULTS 
4.1. Mean He II photionization rate 

We now apply our semi-analytic model to the observed 
He II opacity found in the right panel of Figure [3) For 
each redshift bin, the photoionization rate is calculated 
by iteratively solving equation ([S]), i.e. varying (r) (the 
mean photoionization rate for a fully-ionized IGM) until 
F matches the measurements. As noted in ij'2.21 a mean 
free path (eq. fTU|) is needed to determine /(F) for the 
fluctuating background, but the uniform case has no such 
requirement except insofar as it affects (r). 

First, to provide some intuition, we fix k to be the same 
for both cases (here, K — 0.291) and plot the resulting 
effective optical depth in Figure [4J given the emissivity 
and mean free path from ^2.31 The figure shows the red- 
shift evolution of this optical depth for uniform (solid) 
and fluctuating (dashed) radiation backgrounds. Inter- 
estingly, r e ff is significantly smaller for a uniform back- 
ground, especially at higher redshifts (probably due to 
the shrinking attenuation length assumed in the fluctuat- 
ing model). This is because most points in the IGM have 
T < (r) for the fluctuating background, so most of the 
IGM has a higher opacity; the "proximity zones" around 
each quasar are not sufficient to compensate for this ef- 
fect. The measured opacities are included in the figure 
for reference, showing that the shapes of both (normal- 
ized) models are consistent with observations for z < 2.8 
but deviate significantly at higher redshifts. 

In practice, r e ff is the measured quantity, from which 
we try to infer T. We find that, for the fiducial at- 
tenuation lengths, including a realistic fluctuating back- 
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Fig. 4. — The He II effective optical depth, assuming the HRH07 
QLF and eq. dlOI . as a function of redshift. The normalization 
for all curves is constant, n = 0.291. The solid (dotted) curve is 
based on a uniform radiation background, with 7 = 1.0 (1.6). The 
dashed curve represents a fluctuating background with an isother- 
mal temperature-density relation. The average measured opacities 
are shown for reference. 

ground increases the required (r) by about a factor of 
two - a nontrivial effect that is important for reconcil- 
ing quasar observations with the forest. The magnitude 
of the required adjust ment is comparable to that found 
bv lBolton et al.l (|2006f) . who used numerical simulations. 
However, as we have emphasized above, we cannot use 
our model to estimate the absolute value of the ionizing 
background because the FGPA does not fully describe 
the Lya forest; instead we need a renormalization factor 
k. For the remainder of this paper, we therefore fix the 
photoionization rate at the z = 2.45 HRH07 value (to- 
gether with our fiducial Rq), as described in §2.31 This 
procedure gives k = 0.457 and 0.291 for the uniform and 
fluctuating background, respectively. We strongly cau- 
tion the reader that the remainder of our quoted results 
will therefore mask the overall amplitude disparity be- 
tween these two cases, and we will focus on the redshift 
evolution of T instead. 

Figure [2] displays the mean photoionization rate for a 
uniform (triangles) and fluctuating (squares) UV back- 
ground (normalized to the fiducial point). The curves 
represent the photoionization rate inferred from vari- 
ous quasar luminosity functions, as described in H2.'3{ 
Since the z = 2.45 reference point was chosen arbitrarily, 
the overall amplitude should not be considered reliable, 
which is emphasized by the spread in the QLF curves. 
The difference between the fluctuating and uniform UV 
background is small and certainly within the uncertain- 
ties. The effect of including fluctuations on the photoion- 
ization rate is not straightforward. Generally, at lower z, 
the fluctuating T is smaller than the uniform result, and 
the opposite is true for higher redshift. The redshift of 
the crossover between the two behaviors depends on the 
amplitude of the measured transmission ratio; a higher 
F decreases this crossover redshift. 

We therefore attribute this effect to changing the char- 
acteristic overdensity of regions with high transmission, 
which is larger at lower redshifts because of the Uni- 



verse's expansion. In this regime, where the underlying 
density field itself has a relatively broad distribution, the 
fluctuating ionizing field makes less of a difference to the 
required (r). Remember, however, that these relatively 
small changes are always swamped by the differing k's, 
and a fluctuating background always requires a larger (T) 
than the uniform case. 

For z < 2.7, the normalized points lie near the HRH07 
curve. The averaged opacity at z = 2.25, which is signif- 
icantly lower, relies on a single line of sight, HE 2347- 
4342, and differs significantly from the trend seen in 
Figure [H As expected, the inferred T fluctuates con- 
siderably over this redshift range. In part, these vari- 
ations are due to the limited amount of data, both in 
the number and redshift coverage of usable quasar sight- 
lines. But the UV background fluctuates considerably, 
especially during reionization (see the f(T) discussion in 
i )2.2|) . For z > 2.8, the calculated photoionization rate 
consistently undershoots the model prediction, possibly 
indicating the end of helium reionization around that 
time: there is much more He II than can be accommo- 
dated by a smoothly varying emissivity or attenuation 
length. 

4.2. Evolution of the attenuation length 

Because the measured quasar emissivity evolves 
smoothly with redshift, the most natural interpretation 
of this discontinuity is in terms of the attenuation length, 
which intuitively evolves rapidly at the end of reioniza- 
tion when He III regions merge together and sharply 
increase the horizon to which ionizing sources are visi- 
ble. Following the prescription for the UV background 
in H2.31 we calculate the mean free path Rq, given the 
HRH07 QLF and the T_i4 from the previous section. 
This procedure amounts to varying the solid curve, via 
Ro, to match the points in Figure^ The redshift evo- 
lution of the attenuation length for uniform (triangles) 
and fluctuating (squares) radiation backgrounds is plot- 
ted in Figure [Bj with equation (JTDJ) as a reference. The 
normalization remains the same as the previous section, 
i.e. z — 2.45 is the fiducial point. 

Similarly to the inferred photoionization rate, the 
points vary about the reference curve for z < 2.7 and 
depart from it for z > 2.8. The uncertainties, which 
are shown only for the uniform UV background (but are 
comparable for the other case), are again quite large. In- 
corporating fluctuations reduces the severity of, but does 
not eliminate, the jump in the evolution of the attenua- 
tion length. Once again the results lie consistently below 
the curve for higher redshifts. From the viewpoint of T 
or Rq, there appears to be a systematic change in be- 
havior above z « 2.8. The marked decrease in Rq that 
is required, by at least a factor of two from the fidu- 
cial model, indicates an important change in the state of 
the IGM. However, as we have described above, a single 
attenuation length is no longer appropriate during reion- 
ization, so in Sj5]we will turn to models of inhomogencous 
reionization. 

4.3. The IGM Temperature- Density Relation 

As discussed in H2.21 the temperature-density relation 
of the IGM is a complicated question and an important 
component of the semi-analytic model. For the majority 
of the paper, we assume an isothermal model, i.e. 7 = 1 





Fig. 5. — Evolution of the helium-ionizing attenuation length 
with redshift. The calculated mean free path generally increases 
with a discontinuity around z = 2.8. The fluctuating (squares) 
background appears to smooth the evolution as compared to the 
uniform (triangles) background. The results are matched to the 
I Bolton ct al. (2006) attenuation length (solid curve) at z = 2.45. 
The estimated uncertainties are shown only for the uniform case 
but are similar for both. 



Fig. 6. — Comparison of the inferred photoionization rate (in 
units of 10 — 14 s _1 ) for two temperature-density relations, given 
a homogeneous radiation background. The triangular points are 
derived using an isothermal model. The square points assume a 
steeper temperature-density relation, 7 = 1.6. The points have a 
small redshift offset for illustrative purposes. The photoionization 
rate computed from the HRH07 quasar luminosity function is plot- 
ted for reference. Both models are normalized to the observations 
at z = 2.45. 



in T = To A 7 . In reality, the temperature may depend 
on the density of the IGM. A further complication arises 
during (and shortly after) helium reionization when the 
IGM is inhomogene ously reheated and subsequently re- 
laxes to a power law dGleser et al.ll2005l : lFurlanetto fc Obi 
I2008bl iMcQuinn et al.ll2009f) . 

To partially address the former issue, we repeat the cal- 
culation of r for a homogeneous radiation background, 
but now with 7 = 1.6, shown in Figure [Sj The difference 
between the two cases is < 30%. Overall, the steeper 
temperature-density relation slightly smoothes the jump 
in r_i4 near redshift z w 2.8 but is insufficient to fully 
explain the observed disconitnuity. As expected, the nor- 
malization differs between the two equations of state: 
k 7=1 .o = 0.457 and re 7 =i. 6 = 0.277 (see Fig. gj). In 
other words, a model with a higher 7 requires a higher 
(r) to achieve the same optical depth. This is because 
a steeper temperature-density relation makes the low- 
density IGM, which dominates the transmission, colder 
and hence more neutral. Overall, then, the temperature- 
density relation and fluctuating ionizing background lead 
to a systematic uncertainty of nearly a factor of four in 
the mean photoionization rate inferred from the He II 
forest. 

5. MODELS FOR HE II REIONIZATION 

We have seen that the observations appear to require 
a genuine discontinuity in the properties of the IGM at 
z « 2.8, although large statistical errors stemming from 
the small number of lines of sight prevent any strong con- 
clusions. This change is often attributed to reionization; 
here we investigate this claim quantitatively with several 
"toy" models for the evolution of the helium ionized frac- 
tion XHoiii- This fraction determines f(T) as described 
in £12.21 Assuming the HRH07 QLF and the attenuation 



length given by equation (jTQI) , 3 we calculate the effective 
optical depth via equation ©. The fiducial point for 
normalizing k remains at z = 2.45. 

Figure [7] shows the effective optical depth r e g and ion- 
ized fraction ^Hein as a function of redshift for five reion- 
ization models. Each scenario is characterized by the 
redshift, zh c , at which XHciii reaches 1.0. From left 
to right in the figure, the curves correspond to zh c = 
2.4,2.5,2.7,3.1 and z Hc > 3.8, i.e. post-reionization for 
the entire redshift range in question. The rate of ioniza- 
tion varies slightly between the models. 4 The measured 
opacities are plotted for reference, including t he recently 
discovered S PSS J2346-0016 at z = 3.45 (| Zheng et all 
l2004al l2008h . Note that we do not estimate any cosmic 
variance uncertainty. 

The effective optical depth evolves smoothly in the 
post-reionization regime, which seems compatible with 
the data below z w 2.8. The zh c = 2.7 model seems to 
fit the data best, but a range of values for the reionization 
redshift would be consistent with the existing data. Note 
that models with zjjc > 3 seem to evolve too smoothly; 
however, none of the curves displays enough of a discon- 
tinuity to match the data completely. Obviously, more 
lines of sight at z > 3 are needed to reduce the wide 
cosmic variance. 

6. DISCUSSION 

3 Again, we note that during helium reionization the atten- 
uation length should be evaluated with reference to individual 
quasars; in that case, i t does not evolve strongly with redshift 
lIFurlan ctto & Oh 2008a|). This consideration will only strengthen 
our conclusions. 

4 The unevenness in the curves arises because generating 
our Monte Carlo distributions is relatively expensive computa- 
tionally, so we only generated a limited number at ShcIII = 
(0.3, 0.5, 0.75, 0.9, and 1.0) for each redshift. 
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Fig. 7. — The effective helium optical depth and He III frac- 
tion for five toy reionization models. From left to right, the curves 
correspond to helium fully-ionized by zpjc = 2.5,2.7,2.9,3.1 and 
post-reionization (or zjjc > 3.8). The He III fraction evolution 
is varied slightly between models. The measured opacities are in- 
cluded for referen ce, including the lo wer limit at z = 3.45 from 
SDSS J2346-0016 llZheng et alj|2004al) . 

We have applied a semi-analyic model to the inter- 
pretation of the He II Lya forest, one of the few direct 
observational probes of the epoch of helium reionization. 
Using simple assumptions about the IGM, the ionization 
background, and our empirical knowledge of quasars, we 
have inferred the evolution of the helium phoionization 
rate and the attenuation length from the He II effective 
optical depth. We averaged the opacity measurements 
over 5 sightlines, which show an overall decrease in r e g 
with decreasing redshift and a sharp jump at z « 2.8 
with the alternating low and high absorption at higher 
rcshifts. After proper normalization, our model provides 
good agreement to the lower redshift data, but - assum- 
ing smooth evolution in the quasar emissivity and at- 
tenuation length - consistently overpredicts (r) above 
z « 2.8. Although the uncertainties are large, these re- 
sults suggest a rapid change in the IGM around that 
time. 

Our semi-analytic model is based the quasar luminos- 
ity function and the helium-ionizing photon attenuation 
length, which are determined empirically. The uncer- 
tainty in these quantities significantly affects our results. 
In particular, the plausible range of the mean EUV spec- 
tral index, 1.6 < a < 1.8, shifts the amplitude of our 
results by a factor of about two. Our model takes only 
the mean value and does not account for the variation 
in a from different quasars. Furthermore, our treatment 
of the attenuation length ignores any frequency depen- 
dence. However, these factors likely only affect the am- 
plitude, not the evolution, of the inferred photoionization 
rate, especially the the jump at z rs 2.8. A steeper red- 
shift evolution for the attenuation length would decrease 
the severity of the feature around z ~ 2.8, but a simple 
power law cannot eliminate it, given our method. 

In calculating T and Rq, we compared uniform and 
fluctuating backgrounds. Although he lium reionization 
is tho ught to be inhomogeneous (see iFurlanetto fc Or] 
l2008al) . the assumption of a uniform background has 



been common. We found that the uniform case pro- 
duces an effective optical depth approximately a factor 
of two smaller for a fixed (r) . Thus, properly incorporat- 
ing the fluctuating background is crucial for interpreting 
the He II forest in terms of the ionizing sources. Further- 
more, we find that the inclusion of background variations 
slightly smoothes, but does not remove, the jump in the 
attenuation length at z « 2.8. A clear change in the IGM 
does appear to occur around this redshift. 

The discontinuous behavior in T and Rq led us to in- 
clude helium reionization in our model through the dis- 
tribution /(r). During reionization, the ionized helium 
fraction determines this distribution, and we studied sev- 
eral toy models for the redshift evolution of xhciii- These 
models suggest zh c ~ 2.7 as the best fit to the data, 
but the statistical uncertainties in T er r are large. We do 
not account for cosmic variance and only consider the 
mean effective optical depth. Our method also makes 
assumptions that are not valid during reionization, e.g. 
a power-law temperature-density relation, but this does 
not appear to affect the discontinuity significantly. 

In fact, the most important caveat to our model is 
the use of the fluctuating Gunn-Peterson approximation, 
which is a simplified treatment of the Lya absorption. 
The approach ignores the wings of absorption lines, pecu- 
liar velocities, and line blending; overall, these effects re- 
quire us to add an unknown renormalization factor (of or- 
der ~ 0.3-0.5) when translating from T to optical depth 
and compromises attempts to measure the absolute value 
of (r). One danger is the possible redshift evolution of 
this factor: we have assumed that it does not evolve, but 
in reality the line structure and temperature of the for- 
est do evolve, especially at the end of reionization. More 
detailed numerical simulations that incorporate both the 
baryonic physics of the Lya forest and the large-scale 
inhomogeneities of helium reionization are required to 
explore this fully. 

Interestingly, if our interpretation is correct then it ap- 
pears that helium reionization completes at zjjc ^ 3. 
This places it several hundred million years after the 
epoch suggested by indirect probes of the H I Lya for- 
est. Specifically, some measurements of the temperature 
evolution of the forest show a shar p jump at z ~ 3.2 
and a shift towar d isothermality (jSchave et al.l 120001 : 
iRicotti et al.l l2000f ) , which has been interpreted as evi- 
dence for helium rei onization (|Furlanetto &: Ohl [2008b; 
iMcQuinn et al.l 120091 ). It is not clear whether this time 
lag can be made consistent; it probably depends on the 
details of the line selection in the obs ervations (see the 
discussion in IFurlanetto &: Ohl l2008bh . Moreover, late 
helium reionization would present further difficulties for 
an explanation of the z ~ 3.2 feature in the H I Lya 
forest opacity in t erms of helium reionization (see also 
iBolton et al.l 120 09). Another indirect constraint is con- 
sistent, however, with our picture: reconstruction of the 
ionizing background from optically thin metal systems 
finds an effective optical depth in He II Lya photons 
slightly higher than the di rect measurements, but with a 
similar redshift evolution ([Agafonova et al.ll2005l [2007h . 

The most significant limitation in the data is the rela- 
tively small number of lines of sight, producing large vari- 
ations in the measured transmission, especially at z > 3 
where the cosmic variance is large. Fortun ately, a num- 
ber of new lines of sight have been found (|Zheng et al.l 
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[2008HSvph ers et al.ll2009f ). and the recent installation of 
the Cosmic Origins Spectrograph (COS) on the Hubble 
Space Telescope adds a new instrument to our arsenal. 
Although the nominal wavelength range of COS limits it 
to z > 2.8, this is precisely the most interesting range for 
studying reionization. Our models show that r e ff < 5 so 
long as XHciii ^ 0.3 (see Fig. [7]), so there should be a rel- 
atively wide redshift range with measurable transmission 
- especially when considering the wide variations in the 
ionizing background expec ted during and a fter helium 
reionization (see Fig. [Tl and lFurlanettoll2008l ). 

Finally, these prospects point out one important dif- 
ference between He III and H I reionization: the near- 
uniformity of the ionizing background at the end of H I 
reionization means that very little residual transmission 
can be expected at z > 6 for that event, making the 
Lya forest relatively useless for studying reionization. 
In contrast, the large variance intrinsic to the He II- 
ionizing background produces much stronger fluctuations 
and makes the epoch of reionization itself accessible with 
the He II Lya forest. 

Another interesting difference between helium and hy- 
drogen is the effect of including fluctuations on the pho- 
toionization rate inferred from the Lya forest. We find 
that assuming a uniform ionizing background underesti- 
mates r by up to a factor of two, while during hydrogen 
reioni zation the effect is much smaller - only a few per- 
cent ()Bolton fc Haehneltl 120071 : iMesinger fc Furlanettol 
2009). During and after helium reionization, the fluc- 



tuations are much more pronounced than the hydrogen 
equivalent, leading to a much broader /(J), so that more 
of the Universe lies significantly below the mean. For hy- 
drogen reionization, the distributions are much narrower, 
favoring T near the mean. In addition, after hydrogen 
reionization the density distribution is much wider than 
/(r), so that the latter provides only a small perturba- 
tion; the opposite is true in our case. 

Our gene ral approach is very similar to iFan et aTl 
(2002) and I Fan et all (|2006h . who also interpreted the 
H I Lya forest data at z ~ 6 using a uniform ionizing 
background and the same IGM density model as we have 
(although in their case that model required extrapolation 
to the relevant redshifts). They also found a disconti- 
nuity in the optical depth (at z ~ 6.1), which is often 
taken as evidence for H I reionization. But during this 
earlier epoch, that inference is less clear because of the 
near saturation of the forest and the unknown attenua- 
tion length (whose evolution really determines the overall 
ionizing backgrou nd, but which may evolve rap idly even 
after reionization: iFurlanetto &: Mesingerll2009( ). Never- 
theless, we hope that understanding this discontinuity in 
the He II forest properties will shed light on the problem 
of hydrogen reionization. 

We thank J. S. Bolton, J. M. Shull, and J. Tumlinson 
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